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Abstract 

At HERA, the large flux of almost real photons accompanying the electron beam 
leads to the copious photoproduction of jets. Regions of small momentum fractions 
^ I X of the incoming particles are explored, where the density of partons is high. As 
^ I a result, the probability for more than one hard partonic scattering occurring in a 
single 7p collision could become significant. It is well known that this effect modifies 
O . the contribution of jets (minijets) to the total cross section. We discuss the latest 
^ . HERA data on the total jp cross section in this context. The possible effects of 
^ ! multiple hard interactions on event shapes and jet cross sections at HERA have been 
CL|| studied using Monte Carlo programs. We review some of the available results, which 
p • in general indicate that the effects of multiple interactions should be significant and 
^ ■ may already be manifest in the existing HERA data. 

> ■ 

•rH , 1. Introduction 

X 

5r , For both protons and photons, QCD predicts a rapid increase in parton densities at 
low X, where x is the fraction of the beam particle's momentum which participates 
in the 'hard' scattering (interaction). In a naive treatment, this rise can lead to a 
corresponding (but ultimately unphysical) rise with increasing energy of perturbative 
QCD calculations of the jet contribution to the total cross section. However, the large 
number of small x partons contributing to jet production can mean that there is a 
significant probability for more than one hard scatter per 7p interaction. Only the 
resolved part of the photon can undergo multiple interactions, i.e. the direct part is 
taken to be unaffected. The effects of multiple interactions can provide a mechanism 
for taming the rise in the QCD cross section in accord with unitarity, as discussed 
in Sjostrand's talk p|. 

Additionally, there are clearly implications for the hadronic final state. In order to 
study these effects eikonal models have been implemented within the hard process 
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Figure 1: An example of a multiple scattering in a jp collision. 

generation of several Monte Carlo programs P, |, ||, Multiple parton scattering, 
as illustrated in Fig.l, is expected to affect jet rates. The average number of jets 
per event should be increased when partons from secondary hard scatters are of 
sufficiently high px to give jets in their own right. In addition, lower pt secondary 
scatters produce extra transverse energy in the event which contributes to the pedestal 
energy underneath other jets in the event. Thus multiple scattering can influence jet 
cross sections even when no parton from the secondary scatters is itself of a high 
enough px to produce an observable jet. By boosting the transverse energy of jets in 
this way, multiple scattering leads to an increase in jet cross sections for jets above a 
certain Ej^ cut, even though the total cross section is reduced. 

The theory of multiple scattering is still not well understood; therefore phenomenology 
is based on models with several assumptions and unknown parameters. The effects 
of multiple scattering on the total cross section and on jet/event shapes probe these 
unknowns in somewhat different ways. Although they are related, it therefore makes 
sense to consider the two aspects separately. In sect. 2 we outline the framework of 
models of multiple interactions, in sect. 3 we discuss the HERA total cross section 
data for jp quantitatively from the minijet point of view. In sect. 4 we present some 
of the expected effects of multiple interactions on jet cross sections, with reference to 
existing data, and in sect. 5 we give some conclusions. 

2. Modelling multiple interactions 

When speaking of multiple scattering (alias multiple interactions), the basic build- 
ing blocks are the standard 2^2 partonic interaction processes, such as gg gg, 
qq qq, qq gg and qg —>■ qg. An event may contain none, one, two or more 
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such interactions; the novelty of the multiple interactions concept is that the 'two 
or more' occurrence may be the norm at high energies rather than a rare exception. 
The correlations between interactions occurring in the same event may be horribly 
complicated, but the hope is that it is feasible to approximate this by a simple fac- 
torized description in terms of several independent scatterings. This may be made 
plausible by viewing the incoming hadron and hadronlike photon as two pancakes, 
flattened by Lorcntz contraction. When these pass through each other, different 
parts are causally separated, so the probability of an interaction between any pair 
of partons can be assumed independent of what occurs anywhere else. This gives a 
poissonian distribution in the number of interactions, for a fixed impact parameter 
between the two colliding hadrons/photons. In central collisions the average number 
of interactions should be larger than in peripheral ones; it is therefore necessary to 
introduce a model of how partons arc distributed in the pancake. The 'cikonalization' 
procedure, so central to current descriptions, is a combination of these two concepts: 
a poissonian distribution at each fixed impact parameter plus an impact-parameter- 
dependent overlap function. This form may be integrated over the impact parameter 
to give the total cross section, the probability distribution of interactions, and even 
the number of additional interactions underlying a hard interaction. 

Several objections can be raised. For instance, scatterings are assumed to be disjoint, 
i.e. a parton does not undergo more than one scattering. This can be motivated 
by simple counting arguments: if each hadron contains partons, the rate of two 
disjoint 2 — > 2 scatterings is proportional to but that of two scatterings with a 
shared parton, i.e. a 3 ^ 3 process, is only proportional to A^^. A counterexample 
would be the occurrence of 'hot spots' caused by the cascading of a single parton. 
This is especially relevant since each 2 — > 2 subprocess is embedded in a larger system 
of initial- and final-state radiation. Colour correlations between interactions are not 
well understood, and yet this is of large importance for the structure of the final 
event. In summary, current cikonalization approaches can only be viewed as sensible 
first approximations, and therefore the existence of several different models is an 
advantage. 

Maybe the most ambitious approach is the DTU (dual topological unitarization) one. 
Here a hard pomeron term is given by scatterings with pt > PTmim ^ soft pomeron 
term includes nonperturbative scatterings at small pxi and triple- and loop-pomeron 
graphs are added to incorporate also diffractivc topologies in the description. The 
PTmin scalc Can be changed freely over some range, since the nonperturbative pomeron 
term can be modelled to compensate for the change in the hard interaction rate. This 
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approach is represented by the Phojet program 0. 

In less ambitious approaches, the issue of diffractive events is kept separate, and 
unitarization within the soft sector is simphfied to allow (at most) one soft interaction. 
Here PTmm is a real parameter of the theory, to be determined by a comparison with 
data. Of course the use of a sharp cutoff in pt is only an approximation to what has 
to be a smooth turn on in reality. Examples include the models found in Pythia 
and HERWIG 0, These two are similar in philosophy but differ in the details. For 
simplicity we therefore concentrate on the latter in the following. 

Let us start by considering a jp interaction at some fixed impact parameter, b, and 
centre-of-mass energy, s. In particular, we suppose the 7 to be hadronlike (i.e. re- 
solved 0). We assume that the probability for the photon to interact in such a state 
is Pres and take the p-dominance form, i.e. Pres = 47raem//p ~ 1/300. The mean 
number of jet pairs produced in this resolved-7-p interaction is then 

(n(6, s)) = £partons ® an (1) 

where £partons is the parton luminosity and an is the cross section for a pair of partons 
to produce a pair of jets (i.e. partons with pt > PTmin)- 

The convolution is because the parton cross section depends upon the parton energies. 
More specifically, 

'^•^partons ^^{b^Tl'yi^X'y^Tlpl^Xp^dx^ciXp (2) 

where A{b) is a function which specifies the distribution of partons in impact param- 
eter. It must satisfy 

J ndb'^A{b) = 1 

in order that the parton luminosity integrated over all space is simply the product of 
the parton number densities. Factorizing the b dependence like this is an assumption. 
In particular we do not contemplate QCD effects which would spoil this, e.g. perhaps 
leading to 'hot spots' of partons. Also ni{xi) is the number density of partons in 
hadron i which carry a fraction Xi of the hadron energy. For ease of notation we do 
not distinguish between parton types and have ignored any scale dependence of the 
number densities. For the proton, the number density is none other than the proton 
parton density, i.e. np{xp) = fp{xp). However, since we are dealing with resolved 
photons, the number density is related to the photon parton density by a factor 
of Pres, i-6. ny{x^) = f'y{xj)/Pres- Thus, after performing the convolution, we can 
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write: 

{n{b,s))^^aT{s), (3) 

-'res 

where crff'^{s) is the inclusive cross section for 7p to jets. Restoring the parton indices, 
it is given by 

<^T{s) = [Tdp'r [dx, [\ dx,Y.hh{^..Pl)h/A^vyT) ^^^^^''ff^P^y (4) 



In order to investigate further the structure of events containing multiple interactions 
we need to know the probability distribution for having m (and only m) scatters in a 
given resolved-7-p event, P^. In order to do this we assume that the separate scatters 
are uncorrelated, i.e. they obey poissonian statistics. Thus 

° -'<"'''^>>*-exp(-WM)». (5) 



m! 



This formula is central to the Monte Carlo implementation in HERWIG. 

We can now ask for the total cross section for 7p — > partons with pt > PTmin- 

/oo 
dh^ 5Z Prn 
m=l 

= nPresjdb''[l-eM-Hb,s)))]. (6) 



Since the total inclusive cross section (cx)^'^) counts all jet pairs (even ones which occur 
in the same event) we expect it to be larger than an by a factor equal to the mean 
number of multiple interactions per event (i.e. averaged over impact parameter). This 
is easy to see. Let (n(s)) be the average number of jet pairs produced in resolved-7-p 
events which contain at least one pair of jets, then 



{n{s)) 



Sdh^Y.^^,mPm 

!dh^TTn=lPra 

Idb\n{b,s)) 
J db'^[l — exp(— (n(6, s 
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Note that an must always be less than the total ^yp cross section, whereas cr)^ need 
not be. 



In order to study the details of the hadronic final state in the presence of multiple 
interactions it is most convenient to use a Monte Carlo simulation. The HERWIG 
Monte Carlo generates the required number of hard scatters and the associated initial 
and final state parton showering. The outgoing partons and remnant jets are then 
fragmented to the hadronic final state. Two modifications are made to the simple 
eikonal model in the implementation. Firstly, energy conservation is imposed, i.e. 
after the backward evolution of all the hard scatters in an event, the energy remaining 
in the hadronic remnants must be greater than zero. Secondly, if during the backward 
evolution of the first scatter the splitting gg ^ 7 is arrived at before the evolution cut 
off scale, the event is classified as an 'anomalous' event and no multiple interactions 
are allowed. 



3. Total 7p cross section 



The 'minijet' contribution, anis), is clearly only one contribution to the total cross 
section. To compare with experiment we must add the soft cross section cr^^^, which is 
assumed to contain all of the physics below PTmin- Furthermore there is a (small) cross 
section, a^'^, for the production of single jet pairs (with px > PTmin) by unresolved 
(i.e. direct) photons. Often it is assumed that PTmin = PTmin, but this differs between 
models. Of course simply adding together cross sections is only a first guess. 



As an example of multiple interaction approaches to the total jp cross section, we 
take the eikonalized minijet calculations for ctj/(s) of ref. [§, which used the parton 
distribution functions of refs. |T0| for the proton and photon respectively, and add to 
them various choices of parametrisations of a^^^. We adjust prmin and the parameters 

This is similar 



of (T^-^^ accordingly, and compare to the recent data on cxj? 
to what was done in ref . |]13| , except there the old HERA data were used [ll4| , p!5[| . For 
a fuller discussion of this approach and its theoretical uncertainties see refs. p|, M, |13 



7P 



0, p 



The recent results for the total 7p cross section at centre of mass energies of ^ 200 
GeV are as follows. The HI collaboration {{W) = (^s) = 197 GeV) find a?,^ = 
(165 ± 2.3 ± 10.9) /ib and ZEUS {{W) = 180 GeV) obtain = (143±17)/ib. In the 
HI case, the first error is statistical and the second systematic. The ZEUS error has 
the systematic and statistical errors added in quadrature. 
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Figure 2: Energy dependence of the total 7p cross section. See text for the description 
of the curves. 

In the simplest variant, a^^^ is taken to have the form: 

af\s)=A + B/y/s (8) 

where A and B are constants. This ansatz implies that the entire rise in a^? at 
high energies is due to minijets. It has very little theoretical motivation and makes 
rather an ad hoc separation between hard and soft physics, endowing pTmm with great 
physical significance. In fig. 2 we show as the dotted curve the best fit to the data 
that can be achieved with such an ansatz added to the results of ref. for o"h(s), 
calculated with prmin = 2 GeV . We take A = 107.98 /xb and B = 54.34 jjh GeV. This 
leads to an unconvincing description of the the data in the 10-18 GeV range [|T6|, 
which show evidence of a rise with energy. If one attempts to attribute the low energy 
rise entirely to the minijet contribution then one would be forced to choose very low 
values of PTmin, in the 1-1.5 GeV range (see ref.0), which would lead to much too 
high a prediction at HERA, as was found in ref . [|l^ . This is illustrated by the broken 
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curve in fig.2, where we take a constant a^"'' = 114/ib and prmin = 1-5 GeV in anis). 
Again the need to fit the low energy rise leads to a dangerously high cross section at 
HERA energies (no attempt was made to fit the data below ^ 10 GeV: any attempt 
to do so must worsen the fit to the higher energy data). 

Also in fig.2 we show (the solid line) the results of ref.|0] added to the following form 

c soft 

for a J, ■' 

a^"f\s) = As' + B/y/s. (9) 

The choice A = 78.4 /xb, e = 0.058 and B = 117.05 /ib GeV, with p^mm = 3 GeV in 
(Th(s), ensures an excellent fit to the low energy data (i.e. W < 18 GeV) |[T6[] . 



This latter approach clearly provides the best description of the data, but is slightly 
unorthodox in attributing only part (around half) of the rise in cross section to 
minijets. We could argue that this is a rational view of the situation, less extreme 
than attributing all of the rise to a^^^, as in the soft Pomeron approach [Q. We 



emphasize that the minijet question is not a stark choice between the two extremes of 
either minijets providing the entire rise with energy of the total cross section or being 
absent: indeed, jets are undeniably produced in both pp |l^ and 77 |2^, as well as 
7P, reactions. The important issue is whether they make a significant contribution to 
the total cross section at existing energies. In this picture, with pTmin = 3 GeV, they 
do; anis) being around 20 /ib at -^/s = 200 GeV, although this is very sensitive to 
the choice of pTmm- The increase of the soft cross section could still have its origin in 
multiple interactions in the region < PTmin, but then without a simple perturbative 
description. 



4. Hadronic final state 



In this section we illustrate the possible size and nature of the effect of multiple 
interactions on some aspects of jet cross sections and the hadronic final state at 
HERA. Here the experiments typically measure jets with transverse energies as low 
as 6 GeV. To simulate these events with Monte Carlo models, they use values for the 
minimum transverse momentum of a hard scatter of around 2 < pTmin < 3 GeV. In 
addition, cuts on the 7J9 CM energy are made, forcing it to lie typically in a range 
of around 120 GeV < -y/i < 265 GeV. For these choices the mean number of hard 
scatters per event was found to be at least 1.04 (higher for the lower pTmin values), 
taking the MRS pT| proton and GRV [E2] photon parton densities. 



For the events we generate, since the cross section falls rapidly with increasing pt, we 
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Figure 3: The jet profile in rj, for a) x°^^ < 0.75 and b) x"'"^ > 0.75. c) 
x"'"^ distribution for dijets, at HERA energies, with direct contribution. Including 
multiple scattering (solid line) and with multiple scattering turned off (broken line). 



expect that most of the multiple interactions will have pt ~ PTmm and so their main 
effect will be to increase the mean number of observed jets by boosting the underlying 
Et in the event. Events which contain multiple interactions with high enough to 
be observed as jets in their own right are relatively rare, but their observation (they 
appear as pairs of back-to-back jets) would provide striking evidence in support of 
the existence of multiple interactions. 



Jet finding was performed using a cone algorithm with a cone radius R = 1 |]23|. Jets 
have Et > 6 GeV and pseudorapidity —2 < rj^^^ < 2. The additional Et can be 
seen directly in the jet profile, Fig.3(a), where the Et in the jets is plotted against rj 
relative to the jet axis. The pedestal energy is increased with the inclusion of multiple 



interactions. A variable found to be useful at HERA for distinguishing between 



direct and resolved photon events is the 'observable x^\ It is defined by 



and the sum is over the two highest Et jets in the event. It is the fraction of the 
photon's energy manifest in the two jets of highest Et-, and in leading order is exactly 
the fraction of the photon energy which enters the hard scatter, i.e. direct events 
are peaked at x'^" = 1 whilst resolved events have x"^'^ < 1. The x^ cut isolates 
the resolved interactions, and as expected multiple interactions have no effect in the 
direct case. Fig. 3(b). The extra Et enhances the inclusive jet rate around r^-''^* = 1, 
and an increased sensitivity to multiple scattering can be seen in the 77-^^* and Eij^* 
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Figure 4: Jet profile. The solid circles are uncorrected ZEUS data. Uncorrected 
transverse energy fiow seen in the calorimeter around the jet axis {dET/dSr]) where 
6rj = Tj'^'^^^ — 77-'^*, for cells within one radian in of the jet axis, for jets in the range 
5 GeV< E^^^ < 6 GeV and with < r]^^* < 1. 



distributions in dijet events, which are shown in ref.0. An enhancement at high t]^'^^ 
and low is seen. All plots in Fig. 3 have the direct contribution generated using 
HERWIG included, hence the rise at x"^'"^ ~ 0.8 in Fig. 3(c). 

It can be seen from Fig. 3(c) that the effect of multiple scattering is greater in the low 
^obs j-egion. Multiple interactions are more likely to occur here as it is in this region 
that the higher parton densities occur; also the energy conservation constraint is less 
restrictive. 



In Fig. 4 various Monte Carlo models are compared to preliminary ZEUS data ||25|| . 
The jet profile is shown, and the data are not corrected for detector effects. The 
simulated events have been passed through a full simulation of the experiment. Jets 
have 5 GeV < E^^'^ < 6 GeV and pseudorapidities in the range < t]^'^^ < 1. These 
cuts remove the effect of the forward edge of the calorimeter acceptance and allow us 
to study the effect of the different models on the jet profile independently of the Ej^* 
and 77-^^* distribution of the models. The comparison with the various models confirms 
the general conclusion made by HI that introducing multiple interactions can 
improve the agreement between Monte Carlo models and the data. However, there is 
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a large amount of freedom in the models and, in particular, the effect of introducing 
multiple interactions depends strongly on the parton density of the photon (compare 
the histograms for the multiple interaction models using LACl and GRV). For further 
comparisons with HERA data, the reader is referred to the talk of Steve Maxfield P? 
and references therein. 



5. Conclusions 



Multiple scattering must occur at some CM energy, in order that unitarity is not 
violated: the relevant point here is whether it is present at HERA. Clearly answering 
this question using total cross section data is going to be very difficult, given the 
uncertainties in the soft physics. This point was made very clear by our discussion of 
sect. 3, and in ref. 0. 

However, jet cross sections seem to be a more hopeful place to look. In sect. 4, the 
effect on the hadronic final state of multiple parton scattering in 7p interactions has 
been simulated by interfacing an eikonal model of multiple parton interactions with 
HERWIG. Models indicate that the effect of multiple scattering is significant at HERA 
energies. For reasonable experimental cuts, the inclusion of multiple scattering leads 
to significant changes in inclusive and dijet cross sections which should be understood 
before attempting to unfold to parton distribution functions. 
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